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Abstract. We consider a lattice of weakly interacting quantum Markov pro¬ 
cesses. Without interaction, the dynamics at each site is relaxing exponentially 
to a unique stationary state. With interaction, we show that there remains 
a unique stationary state in the thermodynamic limit, i.e. absence of phase 
coexistence, and the relaxation towards it is exponentially fast for local ob¬ 
servables. We do not assume that the quantum Markov process is reversible 
(detailed balance) w.r.t. a local Hamiltonian. 


1. Introduction 

Stochastic Markov dynamics is a crucial modeling device for (non-equilibrium) 
statistical mechanics, as well as a technical tool for simulation purposes, i.e. Monte 
Carlo algorithms. In the past years, Markov evolutions have been also considered 
for many-body quantum systems, see e.g. [J [51 |T51 [151 [13 IHl [HI IMl [201 I22j- 
A lot remains to be clarified about their physical origin and justification, but in 
the present paper we take a top-down approach and we assume that a Markovian 
dynamics is given. Its generator (Lindblad operator) is a sum of quasi-local terms. 
We allow this dynamics to be non-equilibrium in that it does not need to satisfy a 
detailed balance condition, i.e. we do not ask that it is self-adjoint w.r.t. to a KMS 
state (Gibbs state). We are motivated by the following question: When does the 
Markov dynamics allow for phase coexistence, i.e. multiple stationary states in the 
thermodynamic limit? A property that appears often together with uniqueness, i.e. 
no phase coexistence, is (uniform) exponential relaxation to the stationary state, 
meaning that expected values of local observables approach their stationary values 
exponentially fast independent of the volume and boundary conditions. This notion 
and the related property of ‘rapid mixing’, which was important in recent work 
[5l[6l[l3! on similar questions as ours, is made more precise in the main part of the 
text. 

As such, the question of phase uniqueness is also interesting for classical sys¬ 
tems. One could conjecture that for any detailed balance dynamics in the unique¬ 
ness regime (e.g. at high enough temperature for Glauber type dynamics, etc.) the 
uniqueness is stable against small non-detailed balance perturbation. To our knowl¬ 
edge, this has yet been proven in any generality. Nevertheless, there are numerous 
results on uniqueness and uniform exponential relaxation (and log-Sobolev inequal¬ 
ities), in particular in the regime of ‘complete analyticity’. We refer to works like 
[II1[T1[T3[22] and references therein. 

On the quantum side, uniqueness and exponential approach have been estab¬ 
lished in [211122] for small perturbations of product dynamics with specific gen¬ 
erators of the form Ex — 1, where Ex maps arbitrary operators into operators 
that act as (multiple of) the identity within the spatial set X. In [2], exponential 
relaxation is proven starting from a log-Sobolev inequality, and in [28) . exponential 
relaxation was obtained for a class of cellular automata. 
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Our work is concerned with general perturbations of product dynamics (a clas¬ 
sical analogue was treated in m)- We show that the perturbed dynamics remains 
exponentially ergodic with a unique (infinite-volume) stationary state, see Theorem 
|3.1[ Additionally, if the dynamics depends analytically on some parameter, then so 
does the stationary state. Boundary conditions do not affect the thermodynamic 
limit of the dynamics and stationary state, which is the content of Corollary |3.4| 

Acknowledgements. It is a pleasure to thank Christian Maes for interesting dis¬ 
cussions. WDR and MS are thankful for financial support from the DFG (German 
Research Foundation). NC is supported in part by ISF grant No. 915/12. 

2. Definition of the Dynamics 

Quantum Spin System. We consider quantum spin systems on the the d-dimensional 
lattice d > 1. The lattice distance between two sites x,y £ is denoted by 
d{x,y) and we use the same notation for the minimal distance between two sets of 
sites. For every site x £ let Hx be a copy of a finite dimensional Hilbert space 
TL. The operators Ax ■= B{'Hx) equipped with operator norm || • || constitute the 
algebra of single-site observables. We use the symbol d to indicate finite subsets. 
The physical Hilbert space and the algebra of observables for our quantum spin 
system in a volume A (g Z'^ are then given by the tensor products 

■= (2-1) 

For S C A each A £ As can be embedded canonically in Aa as the local observable 
A® 11, which, we say, has support in S. In the same way also each (super-)operator 
O £ B{As) will be identified with a local operator on B{Aa) — B{As) d H(Aa\s) 
without explicitly mentioning it. The quasi-local algebra A is the C* algebra defined 
as the operator norm closure of the inductive limit Ua-^a with subsets A (g Z^^ 
ordered by inclusion. A state is a normalized positive functional on a C* algebra. 
The operator norm on B{Aa) is again denoted by || • ||. Note that the norm of a 
local operator O £ B{Aa) with support in S' C A may increase with the ambient 
volume A (g Z'^, which motivates the definition of the completely bounded norm 
||0||cb := supa|| 0|| on the space of local operator^ 

Unperturbed Product Dynamics. For each x £ Z'^ let G{x) £ B(Ax) be the generator 
of a semigroup t > 0, of completely positive and identity preserving operators 

on Ax, i.e. a quantum Markov semigroup (QMS) in the Heisenberg picture. 

Our main assumption is that 0 is a simple eigenvalue of each G(x) with cor¬ 
responding rank one projection Qx and that there exist constants M, g > 0 such 
that 

sup||ee(")‘(]l-Q,)||^j^<Me-s* (2.2) 

xei.‘‘ 

for all times t > 0. 

As an immediate consequence, there is a (unique) state gg(x) on Ax such that 

hm a(e*^(")A) =pg(,)(A) (2.3) 

for all initial states a and observables A £ Ax- In such a situation, we shall say that 
the semigroup is relaxing to a unique stationary state. For a volume A (g Z'^, 

the non-interacting dynamics on A is generated by the sum of local generators 

gA:=J^G(x). (2.4) 

xGA 

^Since however ||0||cb ^ clim('Hs)||0||, where the right hand side means the norm in 
which is shown for example in the textbook |24| (Proposition 8.11), we use the the completely 
bounded norm mostly to facilitate volume-independent notation. 
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The corresponding semigroup is relaxing to the unique stationary state 

Qq,k ^xgaSG(x)- (2-5) 

Remark 2.1. Every QMS generator Q on Ax can be written in the Lindblad form 
g{A) = i[H,A]+Y.Kl,AK^-\{K*^K^,A}, A&Ax, ( 2 . 6 ) 

ra 

with a self-adjoint H G Ax and a finite family of (Kraus) operators Let 

K. be the algebra generated by this family. If all G^x) are merely local copies of a 
generator G as above, then the irreducibility property Kip = Hx, for all ip G Hx, is 
a convenient sufficient condition for our main assumption. See e.g. the introductory 
part of [8] for more details. 

The Perturbation. The perturbation is defined by a family of local QMS generators 
{V(r)}, r (E V(r) e B{At). We assume that V(r) = 0 whenever T is not a 
connected subset of lA or if T is the empty set and that the norm of V(r) decays 
exponentially in the size of T: 

|||V|||i:= sup^eiri/'||V(r)|Ub<e, (2.7) 

for some decay length I > 0. The constant e may be viewed as the second natural 
energy scale in our setup besides the decay rate g. The above condition implies that 
the sum of all interactions affecting a single site x is bounded in norm by 
The proof of our main theorem indeed does not require that each local term V(r) 
generates a QMS, but only that it annihilates the identity, and that globally 

Va := ^ V(r) (2.8) 

rcA 

is a QMS generator for every volume A. 

Boundary Conditions. Boundary conditions are implemented by a family of QMS 
generators {>V(r)}, T (g Z'^, yV’(r) S B{Ar)- Again it is assumed that W(r) = 0 
whenever T is not a connected subset of Z'^ or if T is the empty set and that it 
satishes |||W|||b < oo for some (not necessarily small) decay length & > 0. 


3. Results 


3.1. Open Boundary Conditions. The weakly interacting QMS in the volume 
A {£ Z^* and open boundary conditions is generated by 

Ca:=Ga + Va. (3.1) 

Let K A lA stand for the thermodynamic limit along any sequence of finite volumes 
A such that every A (g Z'^ is a subset of almost all of them. 


Theorem 3.1. Assume 1/1 > log 2 and define another decay length 

1 


(3.2) 


1/1 — log 2 

Assume further that e < g and choose a decay rate g' > 0 such that e < g — g'. We 
also abbreviate 

K= , (3.3) 

9-9 - ^ 

then the following holds for some constant C > 0 that may only depend on M 
defined in equation (2.2); 
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(1) there is a unique strongly eontinuous quantum Markov semigroup St on the 
quasi-local algebra A, such that 

||e‘^^(A)-5t(A)|| < (i^:-(3.4) 

for all t > 0, A and local observables A with support in X C A. 

(2) St is relaxing to a unique stationary state g^d exponentially according to 

||5t(A) - Q^d{A)t\\ < Xe-®'‘Cl^l||A||. (3.5) 

(3) is relaxing to a unique stationary state qa, for every finite volume A, 
and furthermore 

Qzd = lim QA (3.6) 

AyZ'i 

is the unique weak* thermodynamic limit. 

(4) Truncated correlations decay exponentially in space for the stationary state: 

\g^diAB) - gMA)g^diB)\ < (K - P||||B|| (3.7) 

for all pairs of local observables A, B with support in X,Y d , XnY = 0. 


The theorem is proven below by a basic perturbative expansion of the dynamics 
exploiting the uniform local relaxation of the product dynamics for the temporal 
bound and the exponential decay of the interaction for the spatial bound. 


Remark 3.2. The definitions of the time scale 1/e and spatial scale I are obviously 
not independent. The situation is more transparent for finite range interactions. 
In that case, the constraint 1/1 > log2 is no longer meaningful, since |||V|||i is finite 
for every 1. Our result then can be roughly summarized in the following expected 
form: one can take any g' < g and any V provided that the interaction satisfies 

ll|V|||o<(g-5>-'''/'' (3.8) 

where R is the range of the interaction. 


Remark 3.3. If the generators G{x) and V(r) depend analytically on some param¬ 
eters without violating the assumptions throughout an (open complex) parameter 
domain, then the expansions used in the proof converge uniformly. As a conse¬ 
quence, also the dynamics and stationary state, i.e. ^((A) and gj^d{A)^ for local 
observables A, depend analytically on these parameters. More precisely, this holds 
if the operators G{x) and V(r) uniformly satisfy (2.2) and and as long as they 
annihilate the identity 1. For a non-trivial analytic dependence on the parameters, 
note that and gA cannot be positive for the whole domain. Furthermore, the 
theorem’s decay estimates remain valid throughout the parameter domain, which 
we however can only prove for decay lengths I' that are larger than (3.2 1. The proof 
of the theorem is given for this more general setting and shows that I/(ai/) — 02 , 
for some 01,02 > I, is a lower bound on l/l'. 


Let us also comment on the relation with some interesting recent work on ex¬ 
tended quantum Markov dynamics. In [5], the authors do not only study pertur¬ 
bations of products as we do, but more generally of (locally generated) dynamics 
e*-^^ which satisfy ‘rapid mixing’. In case of a unique stationary state gA for each 
volume A, this condition means that 

||ei£A(^) _ £,a(A) 11|| < CdAI) e"®‘ (3.9) 

holds for all observables A G Aa and with a function CdAI) growing only poly- 
nomially in the volume. As a main result, it is shown in [5] that rapidly mixing 
systems are stable according to 

||e*^^(A) -e*^^(A)]l|| < C'(|A|)(e-h 5(d(A, A))) 


(3.10) 
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for all local observables with support X C A, where is the perturbed generator 
and b{d) > 0 a decaying boundary correction. On the way to this result they also 
prove that rapid mixing implies the existence of the thermodynamic limit of sta¬ 
tionary states. Note however that it is unclear and an interesting question whether 
perturbations of product dynamics as studied here satisfy rapid mixing. It does 
not follow from our result as we only obtain estimates with prefactors that grow 
exponentially in the support of local observables. This question is interesting also 
because it has been established that rapid mixing furthermore implies exponential 
decay of correlations and an area law (for the ‘mutual information’) for the station¬ 
ary state, see 011 US]- 

3.2. Other Boundary Conditions. The result on the existence of the infinite 
volume dynamics St is already contained in the work of Nachtergaele et al. 
applying Lieb-Robinson propagation bounds m for QMS; see also uni for similar 
results in the context of (classical) Markov Processes. In fact, we use their result 
to prove, as a corollary of the theorem, independence from boundary conditions 
for the thermodynamic limit of stationary states. For that purpose we define the 
following QMS generator in volume A' d with boundary conditions outside a 
bulk volume A C A', 

Ca,a’-=J^a'+ Y. ^(r). (3.11) 

rcA'\A 

Corollary 3.4. Under the conditions of the theorem, there are constants C,l > 0, 
such that the following holds for all A C A' d let a be a stationary state for 
e‘'^A-A', then 

|ct(A) - gz4A)\ < (3.12) 

for all local observables A with support in X G A, where the state g^d is defined as 
in the theorem. 

This result may be viewed as a dissipative version of Yarotsky’s result on ground 
states for weak perturbations of gapped non-interacting systems [32]. Here the 
QMS was defined by sums of local generators, and by definition each of them has a 
non-empty kernel containing the identity. This fact simplifies the problem consider¬ 
ably and is somewhat comparable with frustration-freeness when studying ground 
states of local Hamiltonians. Note however we are here, in case of non-equilibrium, 
dealing with generators that may not be self-adjoint. 


3.3. Examples. 


Quantum Ising Model With Dissipation. For this example we consider a chain of 
qubits, i.e. each Hx, a; G Z, is a copy of H = C^. We denote with cr®, i = 1,2,3, 
the usual Pauli matrices, and with |t) and jj.) the spin up and spin down eigenstate 
of cr^ in Dirac’s notation. Additionally, we put 2cr^ = cr^ ± ia^ such that cr^li) = 
|t),cr“|t) = li)- The unperturbed product dynamics is generated by copies of a 
single-site generator 


g(A) = i[ha^, A] + (T^Aa — ^{A,a^a }, A£B{'H.), 


(3.13) 


where h G K is a parameter indicating the strength of a (transverse) field, and 
where the expression in square brackets is the commutator and in curly brackets 
the anti-commutator. The coherent and dissipative part of the generator commute 
with each other. It gives rise to a unique stationary state gQ(x) = |i)(il at each site 
and (2.2) holds /i-uniformly for the projection Q 3 ,( •) = (|| • |j,)]l, for example, with 
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constants M = 4 and g = 1/2. We take a translation invariant perturbation acting 
only on pairs of nearest neighbour sites defined through 

V{{x,x + l)){A) A€A{^^^+i} (3.14) 

where a/, is the local operator acting on site x as and where the coupling strength 
J S K is another parameter. The coherent part of the dynamics is that of the 
transverse field Ising model with Hamiltonian 

H = h^al + J (3.15) 

icGA {ai,£c+l} 


By our theorem, if the coupling J is small enough (e.g. 2J < see also 

Remark 3.2 on hnite range interactions) and for all values of /i, the stationary state 
for the full dynamics generated by /1 a as in (3.1) is unique even in the thermo¬ 
dynamic limit. We chose the particular form of perturbation above to relate to 
the well-known (Hamiltonian) transverse Ising model, which is exactly solvable. 
See also EH, where similar dissipative models are solved exactly. This property 
is however not relevant here: uniqueness holds for any perturbation that is small 
enough according to the theorem. We may, for example, add a small longitudinal 
field oc at each site, then the model is not longer mapped to free fermions by a 
Jordan-Wigner transformation. 

An immediate interesting question would be whether uniqueness persists as J is 
increased or if there is a transition to a phase with multiple stationary states. It 
seems tempting (though a-priori not justified) to make a comparison with ground 
state properties of the Ising model in transverse field, for which indeed a quantum 
phase transition (in the ground state) occurs at a certain critical value of J/h and 
there is phase coexistence above the critical value. This is also very much connected 
to a question that is not even fully resolved in the classical domain: is it possible to 
devise a one-dimensional cellular automaton with ‘non-degenerate’ noise that can 
have phase coexistence, see mm- 


Transport for Weakly Coupled Self-consistent Heat Baths. In our second example 
we start out with a product dynamics describing a quantum spin chain where each 
site a; S Z is in contact with a separate heat bath at temperature T^. Given 
a local Hamiltonian = Hf G this contact is modelled by a QMS with 
unique stationary state Qg(x) = Z~^ exp{—Hx/Tx), the Gibbs state. Its generator 
Q{x) = G{x,Tx) is detailed balance w.r.t. at temperature > 0, i.e. it is 
self-adjoint for the inner product ‘weighted’ with the Gibbs state gg(x)^ 

{AB)x = TT{A*gY)Bgl/:^). (3.16) 

for some 0 < s < 1 and it depends analytically on Tx, see Remark |3.3| The 
perturbation V is a nearest neighbour interaction of the form 


V{{x,x= i[Vx,x+ir], such that [Vx,x+i, HxHx+i] = 0 (3.17) 


supposed to be sufficiently weak, so that the theorem guarantees a unique stationary 
state gA = £'a('/’i, ■ ■ ■ ,Tn), for all volumes A = {1,... ,N} C Z. Using that the 
interaction conserves energy locally, i.e. (3.17), we can write the continuity equation 


0^dtg^{e^^^iHx)) (3.18) 

= £'A(V({a;- l,x}){Hx)) - gA{V{{x,x + l}){Hx+i)) + gj^{g{x){Hx)) 


= jx{Ti ,..., Tat) — jx+i{Ti ,..., Tat) -I- Jx{Ti,..., Tat) 


where jx is the current from site a: — 1 to a: and Jx the current originating from the 
heat bath at site x. The game here is the following: How should one choose the 
temperatures Tj such that the currents J = J 2 ,..., Jn-i from the heat baths in 
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the bulk all vanish (but not necessarily Ji and Jn). Of course, choosing all Tj equal 
is a solution (corresponding to global equilibrium). But there are other solutions, 
e.g. corresponding to a uniform heat current jsc through the chain, 

Ji = jx =-Jn = jsc, x = 2,...,N. (3.19) 

The idea motivating such a procedure is that the coupled heat baths model chaotic 
degrees of freedom along the chain, so that indeed the steady currents from these 
reservoirs should vanish. The non-trivial question is now: How does the profile of Tj 
look like for a nonzero j^c- The problem can be solved for sufficiently small — 
and the result is that the current jsc and local temperature differences T^+i — 
both scale as 1/iV obeying Fourier’s law jsc = k{Tx){Tx+i — T^) for some heat 
conductivity k, see 0. We do not prove this theorem here as it strays too far from 
the main message of this paper. However, some thought shows that the result of 
the present paper allows to establish this, as it shows that the local currents depend 
analytically and exponentially weak on distant temperature parameters. For related 
results on such chains of ‘self-consistent heat baths’ in classical mechanics, see 
e.g. 0. 


4. Proofs 

4.1. Proof of the Theorem. For every volume A (g the QMS can be expressed 
in the form of the norm-convergent Dyson expansion. For any local observable A 
with support in AT C A, X ^ 0, one obtains 


°*'^^(A) = e‘^^(A) + ^ / dti.. .dt„ e(*-‘")®^VAe 




i)®aVa ... 

(4.1) 

Given a subset if C A we define the projection operator 

Ph{E) := ~ Qx)) G) {^xeA\EQ^) (4-2) 

on Ma- Since the one-dimensional range of each Qx is spanned by the identity we 
find that 

vir)PA{E) = o if rn£; = 0 ( 4 . 3 ) 

and in particular V(r)PA(A) = 0 for every F C A. The locality of the expansion 
and the local relaxation assumption (2.2) become accessible by writing out the sum 


of local terms in Va and inserting the identity in the form 

11 = 5] Pa{E) 

EcA 

multiple times in the Dyson expansion. It can be written as 

ft 

, tn. , / 


e*^A(yi) f dsAI*(s,r„,F;„,A), 

T' TP *-'0 


(4.4) 


(4.5) 


where we abbreviated = (ri,...,r„), En = {Ei,..., En+i), with non-empty 
Fi, Pi C A for all f = 1,..., n, and E^+i C A also allowed to be the empty set. 
The above sum runs over all these n-tupels, n > 1, of subsets of the volume A. 
Furthermore we introduced the operator-valued function 

Ht{s, r„, En, A) 

:= f dti ... dt„_i , 

...{V{T,)e*^^^PA{Ei)){A). 


(4.6) 
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For given r„ and En let us define the set 

D = D{T^) := ULifi. (4.7) 

and we state some rather direct observations regarding Ht{s,Tn, En, A): 

(i) Ht vanishes unless 

El C X, Flj+i \ Fi = Fli \ Fi and EinTi^^h, l<i<n (4.8) 
where X is the support of A. 

(ii) Ht does not depend on the volume A as long as X,D C A. 

(iii) Ht does not depend on the time t if En+i = 0. 

(iv) Ht is a local operator whose support is within En+i- In particular it is pro¬ 
portional to the identity if En+i = 0. 

The above expansion, or more precisely the integrand in ( |4.6[ ), can be depicted 
by diagrams that spread into space-time starting from X at time zero and that are 
connected as a consequence of observation (i), see Figure 



Figure 1. (a) Sample diagram for n = 7: the vertical lines indicate 
perturbations V(ri) acting on Fi at time ti, the horizontal zig-zag lines 
indicate the sets within En. 


We continue the proof with two Lemmata; the first one concerns the time inte¬ 
gral of Ht{s, r„, En, A) and then we show A-uniform summability over the spatial 
subsets r„ and En in the second Lemma. % denotes the indicator function. 


Lemma 4.1. For every l',g' gM., g' < g, there there is I" > 0, such that for all 
volumes A d r„ and En within A, all times t > 0, and n>l, 



ds ||i7i(s, Fyj, En, A) 


< exp{-g'tx[En+i 7 ^ ^-g'tix[En+i 


n-nD\} 


(4.9) 
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X 


i-1 


{g-g')\Ei\ 


holds for all times 0 < ti < t 2 < t. Here, one can take I" = V and C = 1 under 
the conditions of the theorem, and l/l” = log(M + 1)(1 + l/l') and C = M + 1 for 
the more general setting with complex parameter dependence introduced in Remark 

EH 


Proof. In the proof we first focus on the more general setting. By observation (iii) 
we can, for every given and E^, restrict the volume to A = X U Z?. Note that, 
for alH > 0 and E C A, 

e*^-PAiE) = (4.10) 

is a product of single-site operators, which are bounded in the completely bounded 
norm by Me“®* or M -|- 1 by assumption (2.21. The integrand in (4.6) can be 
bounded in norm by 


...(V(ri)e‘i^^Px(i^i))(A)|| 
< _|_ x)ir"l g-9(t-s)|r„n£;„+i| 

X II PhiEn)) ■ ■ ■ 

...(V(ri)e‘i^^Px(i^i))(A)|| 
Repeating this step n times gives the upper bound 


(4.11) 


(A/ -h l)l^l+51iiril e-s(‘ 


-s(t-s)|s„+i| 


piiH' 




l|v(r.)||cb 


(4.12) 


where we set to := 0 and as before s := The diagrammatic representation can 
be useful for this exercise. Since |i4| < then get 

f ' ds \\Ht{s,Tn,En,A)\\ exp{g'tx[En+i ^ (ll]+g'tix[En+i = ^ + \D\/l'} (4.13) 

Jti 

<{M + l)l^l+(i+i/«')EJ£;4p|| 


dtr. 


dti... dtn-1 




l|v(roilcb. 




i=l 


By changing the integration variables to := ti—U-i and extending the integration 
domain, we obtain another upper bound for (4.13), 

n poo 


>>- nOO 

(M+l)l^l+(i+i/'')^.l®-l|iA||P[y ds,e-s'*'l^-l|lV(r,)||eb (4.14) 

= (M + 1)I^|||A|| f[ i)^'+'/''^"’-'l!r(C)llcb 


{g-g')\Ei\ 


This finishes the proof of the lemma in the setting of Remark 3.3 


Under the conditions of the theorem, we proceed just as above, but we can 
additionally exploit the fact that HQ^IIcb = 1) since Qx is the t —^ oo limit of a 
(contractive) QMS. □ 
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Lemma 4.2. Recall the definition of K in (3.3) of the theorem. For every I” > 0, 

<{K- l)2l- 


r„,E„:(i) i=l 


e' llV(rj)||cb ^ 1\r,|X| 


{g-g')\E, 
*/ lll'^llli < 9 - 9 ' for l/l = log2 + 1/r. 


(4.15) 


Proof. First recall the constraint (4.8) in observation (i). Fixing Ei and F^ for some 
1 < i < n fully determines a family of at most 2l^’l different compatible subsets 
i5i+i. On the other hand, fixing Ei determines a family of compatible subsets F^. 
Then, for every Ei (g 


||V(ri)||cb 


Ti Ei. 


x^Ei 


( 5 - 5 ') 1-^*1 

2ir.lelFI/0||v(r,)||,b 

{9-9')\Ei\ 


(4.16) 


< JLJE with i/; = iog2 + i/r 

9-9 

where we put primes to indicate that we sum over compatible subsets depending 
on Ei and Ei^Ti respectively. Therefore we obtain the majorizing geometric series 


^z-^eir.l/'"||V(fo) 


cb 




{9-9')\Et 


(4.17) 


^ ||V(ri)||cb el^"l/' llv(r„)||cb 

2^ 2^ 2^ 2^ {g-g')\Ef ■■■ 2^ 2^ 


Ti — \ E\(2X Pi E 2 


Tn Er^ 


{g-g')\En\ 


,,,v,^(jm) 


n—1 


if |||V|||/ < g — g', which proves the Lemma. 


□ 


The above Lemmata together show that the expansion (4.5) converges in norm, 
even if weighted exponentially. 


Eds||i7t(s,r„,L;„,^)|| (4.18) 

X exp{g'tx[-E^n+i 7^ ^+9'tix[En+i = ^]+\D\/l'} 

< (iF- 1)01^111^11 


for all times 0 < <i ,<2 < i and with some constant O > 0 (that only depends on 
M), if e = |||V|||; < g — g', where either 

l// = log2 + l/r or 
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i/i = iog2 + iog(M + i)(i + i/r) ^ ^ 

in the setting of the theorem or of Remark |3.3| respectively. 

In the following, the constant C represents a family of constants and its value 
may change from line to line. Let us describe further consequences of the Lemmata: 
(1) The thermodynamic limit of the dynamics 

5t(y4) : = lim (4.20) 
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ds Ht{s,Fn, En, A) 


r„,E„ 

exists for each local observable A G Ax, X lA, and t > 0. Since the local 
observables are dense in A, the infinite volume quantum Markov semigroup St of 
the theorem can be defined by continuous extension. Let A, A'd lA with A C A'. 
In the difference 


e‘^*(A) - (A) 


,A' 




ds Ht{s, r„, En, A) 


(4.21) 

(4.22) 




all those terms that, diagrammatically speaking, do not reach outside of A cancel 
by observation (ii). Moreover the remaining terms in the sum decay exponentially 
in the distance d{X,A'^) < \D\, 


|e*^A(^) _e*Av(^)|| < (K-i)e-''(-^AA/i'cl^l||A|l. 


(4.23) 


This shows that (3.4) holds by taking the limit A' ^ and that the convergence 
of ( |4.20 1 is uniform for all times. Therefore St is strongly continuous. 

(2) We take g^d to be the state on A defined by 

gzdiA)l= gg,xiA)l+ lim x[En+l=^ f ds Ht{s,T^, E„, A) (4.24) 

f —>-00 Ir, 

for local observables and extend it by continuity to the full algebra of observables. 
We find that 


iSt(A) — gid{A)l 
- e*®^ ( 

x[Eu+i A ' 


(4.25) 


= e-^^- ^A) - pe.x(A)ll 


r„,E„ 


X[En+l = I 


dsiJ4(s,r„,£;„, A) 

ft' 

dsHt'{s,Tn,En,A). 




Due to the fact that the unperturbed dynamics is relaxing exponentially fast, the 
first term is bounded by 


|e*e-(A)-pg.x(A)]l|| = ^ {A^-Px{E)){A) 

li^EcX 

< 2l^l(M +1)1^1 e-3‘||A||. 


(4.26) 


This bound only uses (2.2), which is convenient in view of the remark 3.3 on analytic 
dependence. For products of QMS one can easily improve the above bound to obtain 
a prefactor that grows only linearly in |X|. By (4.18), the other two terms are 
each bounded by 

(AT-l)e"®'‘C'l^l||A||, (4.27) 

so that the theorem’s claim ( |3.5[ ) follows for some appropriate constant C. 

(3) Defining states gx as above, but restricting to finite volumes A, it also follows 
in the same way that e*'^A is relaxing to a unique stationary states exponentially 
fast, 

||e‘^^(A) - pa(A)] 1|| < iFe-9'‘Cl^l||A||. 


( 4 . 28 ) 













12 


N. Crawford, W. De Roeck, and M. Schiitz 


Furthermore, we have 

\qk{A) - < {K-I)e-<^^^A^'^/^'C\^\\\A\\. (4.29) 

It follows that Bid is the (unique) weak* limit of these stationary states. 


(4) Exponential decay of correlations follows from the expansion and (4.18), since 
all terms with |i4| < d{X,Y) cancel in the difference (3.7). 


4.2. Proof of the Corollary. By assumption, the local generators of the bulk and 
boundary dynamics V and W decay exponentially in the diameter of their support 
(assumed here to be connected sets). This infers the existence of an effective ‘prop¬ 
agation speed’ in the system, see [23]: adapted to our purposes we may conclude 
that there exist constants C, /i, u > 0, such that 

||e‘^A,A'(^) (4.30) 

for all local observables A with support in X C A and times t > 0. [To facilitate 
transferring the results from [^ note that, in case of the d-dimensional lattice, 
one may take e“^"(l for the function F^{n) appearing in this reference.] 

Inserting the time evolution appropriately and using the triangle inequality as in 
[Bj we obtain 

|cr(A) - pa(A)| (4-31) 

< |cr(e*^A,A' (y[)^ _ | _|_ |cr(e*^A(^( _ (T^g^(A)]l) I 

< ||e*^A.A'(^) _e‘^A(^)|| ||e*^A(^) _ g,^(A)]i|| 


for arbitrary t > 0. Choosing this time through vt = d{X,A‘^)/2 hence gives 

|ct(A) — pa(A)| < f := max{Z', 2/^} (4.32) 


by the propagation bound (4.30) and exponentially fast relaxation in finite volume 
(4.28). Together with (4.29) we arrive at the Corollary. 
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